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We report the systematic investigations of the oscillation of the Josephson vortex (JV) flow resis- 
tance in Bi2Sr2CaCu208+5 micro-fabricated junctions with various geometries and superconducting 
anisotropy parameters. As the applied magnetic field parallel to the afc-plane is increased, oscillation 
with a period corresponding to a </>o/2 par atomic Josephson junction changes to oscillation with a 
doubled period. This crossover is scaled by both the junction length and the anisotropy parameter, 
indicating that the bulk inductive coupling that favors the triangular JV lattice is replaced with 
the surface deformation energy as the dominant interaction for a JV lattice. These results suggest 
that the in-phase square JV lattice is pronounced at a higher magnetic field in a smaller and more 
anisotropic sample. 



Highly anisotropic high-T c superconductors (HTSC) 
such as Bi 2 Sr 2 CaCu208+,5 (Bi2212) are well described as 
stacks of weakly coupled Josephson junctions (JJ's) that 
consist of CuC>2 superconducting and charge reservoir 
block layers;^ referred to as intrinsic Josephson junc- 
tions (IJJ's). Not only the inductive coupling^ but also 
the capacitive coupling 4 between the stacked Josephson 
junctions are responsible for such peculiar properties of 
IJJ's as multiple branches in the current-voltage (/V") 
characteristics and the longitudinal Josephson plasma 
mode. 5 Rich collective excitation phenomena attributed 
to these couplings among identical stacked junctions^ can 
be adopted for investigating physics of HTSC and also 
promise us potentials for future device applications. 

In the magnetic field parallel to the ab plane, where 
Josephson vortices (JV's) quantized as a unit of 4>o = 
h/2e are induced, various static and dynamic arrange- 
ments of JV's are expected. 7 -^ The variety of phenomena 
is given by the couplings between the adjacent JJ's, which 
make IJJ's a sharp contrast with artificially stacked JJ's. 
The other characteristic features of the JV system in par- 
allel magnetic fields have been found in the IV charac- 
teristics j2il£ the Josephson plasma resonance^ and etc. 

Recently, Ooi et alr^ found a periodic oscillation of the 
JV flow resistance with a period of </>o/2 in Bi2212 IJJs 
with widths of 20-50 /im as a function of the magnetic 
field. They interpreted this is due to the coherent flow 
of the triangular JV lattice. Subsequently, Koshclcv 13 
and MachidaM explained the oscillation by calculating 
the JV flow resistance by taking the surface barrier ef- 
fect into consideration. In this model, since the surface 
barrier introduced by the boundary condition at the junc- 
tion edges only influences the positions of JV's close to 
the surfaces, a bulk property (triangular JV lattice) is 
probed by the surface current. Reducing the junction 
width to ~ /im is expected to make the surface barrier 
effect compete with the bulk inductive coupling between 
the JV arrays in neighboring layers and result in drasti- 
cally different phenomena. In smaller junctions < lO/im, 



the authors found oscillations with the doubled period <po 
at high magnetic fields^ and Hatano et al., reported the 
<f>o oscillation in a sub-micron junction. 1 ^ Motivated by 
our works, Machida 17 and Koshelev 18 proposed advanced 
calculations that the triangular and the square JV lat- 
tices appear alternatively as a function of magnetic field 
and that the square lattice is more favorable for higher 
magnetic field and in samples with smaller widths. Since 
the square JV lattice is a result of in-phase locked os- 
cillation of the Josephson currents of the stacked IJJ's, 
realization of a square JV structure has been considered 
as a central issue for an electromagnetic radiation from 
IJ J. Although recently reported strong terahertz emission 
from IJJ's in the absence of applied magnetic fiel d 19 ' 20 is 
attributed to the size-restricted cavity resonance of the 
rather large mesas (~ 100/im), tuning excited wavelength 
by JV lattice like a field tunable photonic crysta l 11 ' 21 is 
a next goal of the IJJ emitter. 

In this paper, we present the experimental results of 
the JV flow resistance and IV characteristics in Bi2212 
micro-structured single crystals with dimensions of 1—6 
yitm more clearly and systematically than in the previous 
worki 1 ^ In the smaller samples and at higher fields, we 
found periodic oscillation of the JV flow resistance with 
the period of the 4>o instead of </>o/2 as argued by the 
theories. The results clearly indicate that the crossover 
of the two oscillations is well tuned by both the super- 
conducting anisotropy parameter and the junction size, 
which may prove formation of the square JV lattice in 
the static and slowly moving dynamic states. 

Bi2212 single crystals grown with the traveling float- 
ing method were fabricated into bridge-shaped meso- 
scopic junctions with a focused ion beam (FIB) machine 
SMI2050MS (SII NanoTechnology Inc.) .22 The fabrica- 
tion details are described in Ref^. Here, we report the 
results of seven junctions obtained from different slightly 
over-doped crystals as listed in Table [IJ T c 's are in the 
range of 84 - 90 K. External magnetic fields were applied 
by split pair superconducting magnets that generate a 
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TABLE I: List of samples. L, W, and t denote length (± H), 
width (|| H), and thickness (|| c) of junctions, respectively. 
Oscillation starting field H a and crossover field H^/i were 
directly given by measurements. 7 was derived from Hi/ 2 - 
Units for dimensions and magnetic fields are /jm and kOe, 
respectively. 
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horizontal magnetic field up to 80 kOe. Rotating the 
magnetic field around the a6-plane with a high precision 
rotator, c-axis resistance suddenly increases in the vicin- 
ity of the afr-plane.— This resistance is attributed to the 
Josephson vortex flow resistance, and the lock-in state 
without pancake vortices is realized in this angle range. 
We set the field alignment parallel to the a6-plane with 
an accuracy of 0.01 degree. 

Figures HJa) and (b) represent the magnetic field de- 
pendencies of zero-current differential c-axis resistivity 
p c = LW/t ■ dV/dI\i= Q at 60 K in sample B19 and P52, 
respectively. The data were derived from IV character- 
istics measured under external magnetic fields parallel 
to the afe-plane changed step by step. This method for 
evaluation of the JV feature is less ambiguous than the 
method with the ohmic resistance V/I. The p c in P52 
smoothly increases below 7 kOe with increasing field from 
zero and then begins to oscillate with an oscillation cen- 
ter monotonically increasing, as shown in Fig. [ljb). Os- 
cillation period H p is 1.2 kOe, which approximately co- 
incides with adding a half flux quantum 4>q/2 to a JV 
array in the block layer as Hq/2 = 4>$/2Ls — 1.3 kOe 
for P52, where s — 15 Ais the periodicity of the CuC>2 
double layers of Bi2212. This oscillation is quite sim- 
ilar to the previous results*^ The boundary effect for 
the triangular JV lattice is considered to be enhanced 
when the number of JV's in a layer corresponds to ei- 
ther integer (H/Hq = k; full-integer matching) or half 
odd integer(_ff/iJo = k + 1/2; half-integer matching) due 
to matching of the periodicity of the JV lattice with the 
junction width for both two cases. These two matchings 
yield the oscillation with period of Hq/2 and the bottom 
at H/H = k and k + 1/2. 

With increasing magnetic field, the full-integer match- 
ing was diminished above 25 kOe and finally disappears 
above 40 kOe, while the half-integer matching does not 
change to 50 kOe. As a result, H p becomes 2.43 kOe, 
which is twice the value of the H p below 35 kOe and 
close to the Hq value of 2.50 kOe. This crossover from 
(j>o/2 to </>o oscillation was observed in a lower field re- 
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FIG. 1: (Color online)Field dependence of differential c-axis 
resistivity in B19 (a) and P52 (b) at 60 K. Solid symbols rep- 
resent differential resistance dV/dl obtained from IV charac- 
teristics. Broken curves connect maxima and minima points 
in p(H). In the inset of (b), solid curve indicates ohmic resis- 
tances V/I obtained by sweeping magnetic field at / = 1.2/iA. 
Blue vertical grids in (a) and inset of (b) indicate that the 
fields of H/Ho are integers. 

gion in a sample with smaller L. Figure D^a) shows a 
p c — H curve in B19 with L = 1.9 pm. 4>q oscillation 
was certainly observed above 20 kOe, and only a small 
dip at H/Hq = 2 was found as a remnant of the (j)o/2 
oscillation, as shown in the inset of Fig. [lja). 

Since linear IV characteristics were observed even in 
the low current limit as displayed in Fig. [21a), it was 
implied that there is no true critical current in IJJ with 
JV lattices. This makes a remarkable contrast with an 
ordinary vortex system under the c-axis magnetic field 
for HTSC, where a super-linear increase in voltage in IV 
characteristics is regarded as the critical current. In order 
to compare with theoretically-derived field dependencies 
of critical current density J C (H), the current Iq(H) at 
which the induced voltages along the junction are 2 mV 
is plotted as a function of magnetic field at various tem- 
peratures in Fig. 02b). At 75 K below 15 kOe, I (H) 
oscillates with a period of Hq/2 as proposed in Ref. 13 
for J C (H), while Io(H) above 35 kOe shows oscillation 
with the period being Ho as well as p c {H). The behav- 
ior of Iq(H) in the high field region turns to be similar 
to J C {H) in a single Josephson junction (Fraunhofer pat- 
tern) that has local maxima at fields corresponding to 
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FIG. 2: (Color online) (a) IV characteristics at 60 K in V31 
(T c = 84 K) at various magnetic fields from (top) to 60 kOe 
(bottom), (b) Field dependence of Ji, current at V = 2mV 
(the arrow in (a)), at 60, 75, 80, and 82 K from top to bottom. 
Inset shows temperature dependence of anisotropy parameter 
7 obtained by Eq. (JTJ). 



(k + l/2)0o. Such a Fraunhofer oscillation in J C (H) of 
multilayered system like Bi2212 is expected when condi- 
tion L Xj = js is satisfied.— However, Aj is estimated 
to be approximately 0.2 /zm for 7 = 100 — 200, which is 
much smaller than the L of these samples. 

The experimental results demonstrate that the 4>o os- 
cillation is a phenomenon induced by the external mag- 
netic fields, thus it is important to consider competition 
between the bulk and the boundary effects for JV's. At 
relatively low fields where JV's start to form dense trian- 
gular lattices, the JV positions are mostly determined by 
the sharing interaction between vortices in adjacent lay- 
ers due to the inductive coupling between layers. With 
increasing magnetic field, vortices inside the layer be- 
come closer, and interaction between intra-layer JV's is 
pronounced. The JV system is finally dominated by the 
surface deformation energy that lets JV's align in phase 
and turns toward the square lattice in the high field re- 
gion at the cost of sharing energy between JV arrays. 
Since the sharing energy of JV's in a smaller junction 
with less JV's is smaller while the surface deformation 
does not depend on L, the 4>o oscillation is observed at a 
lower field. The single-junction-like behaviors are more 
pronounced in samples with shorter L and at higher H. 13, 



The crossover between the triangular and the square 
JV lattices is quantitatively treated by KoshelevJ^ He 
derived the amplitudes of the resistance oscillation due 
to the half-integer matching and the full-integer match- 
ing, SRi and SR2, as functions of magnetic field. The 
magnetic field at SR2/SR1 — 1/2 is given by the junction 
length L and the anisotropy parameter 7 as: 
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In order to extract generic boundary effect from our re- 
sults in many samples with various oxygen concentra- 
tions, the sample specific parameter 7 should be taken 
from measurable parameters. The comparison between 
the afe-plane resistance before FIB fabrication the c-axis 
resistance^ does not work in this case because oxygen 
would be reduced during FIB fabrication in the vac- 
uum and implanted gallium ions may damage the junc- 
tion. These influences are more serious in junctions with 
widths less than 5 microns. To obtain the anisotropy 
parameter in the exactly same situation as the sample, 
we use the formation of the dense triangular JV lat- 
tice, which is considered necessary for <po/2 oscillation. 
Using a minimum field for triangular dense JV lattice 
H tr = 1.330o/27T7S 2 ,— we obtain a simple relation: 
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where the magnetic fields are given in the unit of Tesla. 
Assuming that H tr corresponds to H s , the bare boundary 
effect for JVs may be extracted from the results. 

Figure [3] shows the values of Hi/2/H^ as a function 
of junction length L. H 1 / 2 /H^ is roughly proportional 
to L. H1/2 and H s were obtained by extracting 8R\ 
and 8R2 from the resistance minima and maxima, as 
shown in the inset of Fig. [3l 8R\ and SR2 are given 
by subtracting R m %n,\, the smooth curves connecting re- 
sistance minima at H/Hq = k + 1/2, and R m in,2 from 
Rmax in p c (H), respectively. Since it is difficult to decide 
the onset of the oscillation, H s was defined as 8R\ = 
by extrapolating R m in,i/ Rmax to 1, so that H s may in- 
clude more ambiguity than U\j2- The slope of the fit- 
ted line is slightly smaller than Eq. ([2]) with a factor of 
0.68, meaning H s = 1.2H tr - Below H tr JV lattices in- 
commensurate along the c axis are expected in the bulk 
as results of the sharing transitions from the triangular 
JV latticed Considering the boundary effect which may 
suppress the sharing transition, the deviation could be 
reasonable. The result that H 1 / 2 can be scaled by 7 sug- 
gests that the crossover between <fio/2 and 4>q oscillations 
is a result of competition between bulk sharing energy 
and the edge boundary effect for JVs. The argument 
that the <f>o oscillation is attributed to the Fiske step at 
higher bias current 26 is not supported because dV/dl at 
1 = shows almost same behavior as ohmic resistance 
V/I, as shown in the inset of Fig. [IJb). Therefore such 
a non-linear effect in voltage should not be considered in 
this current region. 
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with Eq. (JTJ) and plotted it in the inset of Fig. [2jb) as 
a function of temperature. It has been reported that 7 
drastically increases in the vicinity of T c , 27 i 28 It was also 
found that the amplitude of the oscillation of both p c (H) 
and Io(H) is maximized at 60-70 K. This can be under- 
stood that the synchronization of stacked junctions is 
optimized at the temperatures: thermal fluctuation be- 
low 60 K smears the inhomogeneities contained in the 
crystal and helps form regular JV lattice, but thermal 
fluctuation above 70 K smears the matching conditions. 

In summary, we investigated the Josephson vortex lat- 
tice structure in mesoscopic (~ /im) IJJ's with a probe of 
the JV flow resistance along the c axis. The c-axis resis- 
tivity p c oscillates as a function of applied magnetic field 
with a oscillation period changing from (po/2 to 4>q, which 
indicates triangular JV lattice and a possible square JV 
lattice, respectively. Crossover field H 1 / 2 linearly de- 
creases with junction length L. This is attributed to the 
interplay between the boundary effect and the inductive 
coupling between JV arrays in IJJ's. 



FIG. 3: (Color online) L dependence of crossover field i/1/2 
in which 7 is renormalized. Inset: Method to extract H a 
and -Hi/2- Solid and broken curves are those in Fig. [TJb) 
normalized by R m ax (red (dark gray) broken curve). 

The temperature dependence of the crossover field is 
also consistent with the preceding discussions. As tem- 
perature is increased, -ffi/2 considerably shifts to a lower 
field from 60 to 82 K as shown in Fig. [5J We derived 7 
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